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ABSTRACT 

Recent  modern  optimization  algorithms  are  often  very  successful  in  solving  NP-hard  improvement  problems. 
One  of  the  main  disadvantages  of  the  search  algorithm  in  cuckoo  (CS)  is  the  weakness  of  its  ability  as  a  local  search 
algorithm  and  slow  in  its  speed  of  convergence,  which  results  in  a  limited  access  to  the  target.  There  are  two  categories 
that  can  be  solved  by  this  algorithm:  First,  the  unconstrained  functions  are  solved.  Second,  which  is  the  focus  of  this 
article,  is  the  constraint  functions  and  how  to  deal  with  them  using  a  rational  formula  of  the  Barrier  function.  The 
results  show  that  the  CRB  ( Cuckoo  &  Rational  Barrier  Function )  algorithm  has  a  better  convergence  velocity  than  the 
algorithm  (CS)  by  observing  the  separation  of  the  numerical  results  of  the  browser  in  which  the  accuracy  of  the  new 
algorithm  on  some  engineering  functions. 
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1.  INTRODUCTION 

The  main  theme  that  led  to  the  need  to  find  intelligent  search  algorithms  is  to  solve  the  problems  that 
consist  of  multiple  dimensions  as  an  example  of  engineering  problems.  Most  evolutionary  optimization  problems 
in  the  nonlinear  engineering  have  many  constraines.  Metaheuristics  algorithm,  based  on  the  swarm  intelligence,  is 
an  important  part  of  the  global  contemporary  optimization  algorithms.  Good  examples  are  the  genetic  algorithms 
([12],  [10]),  a  model  developed  from  the  Particle  Swarm  Optimization  (PSO)  [14]  and  Harmony  Search  algorithm 
(HS)  [32].  A  number  of  recently  developed  metaheuristics  algorithms  have  been  developed  based  on  the 
intelligence  of  the  swarm  in  nature  ([14],  [7]).  Later,  evolutionary  optimization  algorithms  have  gone  a  long  way  in 
speeding  evolution,  becoming  a  significant  topic  in  improving  the  structure  ([26],  [30]).  In  2009,  both  Yang  and 
Deb  [30]  introduced  the  concept  of  the  Cuckoo  (CS)  search  algorithm,  and  defined  it  as  a  highly  appropriate 
algorithm  for  evolution  and  In  [31]  the  engineering  side  was  introduced  using  the  optimization  method  and  the 
design  of  springs  and  welded  beam  structures.  Day  by  day,  new  algorithms  are  being  developed,  demonstrating 
their  strength  and  effectiveness.  As  an  example,  a  comparison  shown  by  the  cuckoo  algorithm  with  traditional 
algorithms  yields  superior  performance.  The  advantages  of  the  CS  algorithm  include,  first:  enable  it  to  search  the 
world  well  and  second:  the  number  of  parameters  used  in  the  search  is  few.  Gandomi  et  al.  [9].  They  determined 
that  the  best  solutions  they  found  using  CS  were  the  best  solutions  from  solutions  obtained  through  genetic 
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algorithms  and  PSO.  The  CS  algorithm  was  widely  used  in  structural  design  problems  [14],  aerodynamic  improvements 
[17],  and  reliability  improvement  problems  [27].  The  reason  for  the  slow  speed  of  the  CS  algorithm  is  that  it  does  not  use 
gradient  information.  Which  makes  it  difficult  to  apply  this  algorithm  to  complex  engineering  models  where  the  cost  of 
computation  that  exceeds  the  tolerance.  In  [16]  the  researchers  developed  and  evaluated  the  performance  of  three 
algorithms  for  searching  with  Cuckoo  based  on  increasing  the  conversion  parameter  between  local  and  global  random 
paths  as  the  number  of  repetitions  increased.  The  first  algorithm  has  an  increase  in  the  linear  switching  coefficient,  and  the 
second  parameter  is  applied  to  a  steadily  increasing  conversion  as  the  third  parameter  is  used  to  convert  the  increasing 
energy. 

2.  THECUCKOO  ALGORITHM  AND  BARRIER  FUNCTION 
2.1.  The  Cuckoo  Search  Algorithm 

Now  take  a  deep  look  at  the  cuckoo  algorithm  and  explain  its  important  details.  Cuckoo  birds  behave  currish  they 
will  use  other  bird  nests  to  lay  their  eggs.  These  birds  will  be  called  birds  of  the  cuckoo  birds  [19].  Host  birds  usually  know 
when  cuckoo  eggs  are  in  their  nests  either  leave  their  eggs  or  abandon  those  nests.  The  cuckoo  can  adapt  himself  by 
placing  eggs  simulating  host  eggs  in  terms  of  colour  and  pattern.  Giving  them  the  ability  to  reproduce  themselves  and 
reduce  the  number  of  eggs  disposed  of.  As  for  the  ideal  rules  of  the  CS  algorithm: 

•  The  egg  placed  in  each  bird's  nest  is  placed  in  a  randomly  selected  nest  each  time. 

•  Future  generations  will  continue  with  the  best  nests  (good  eggs). 

•  The  egg  developed  by  the  cuckoo  is  inspected  by  host  birds  of  a  fixed  number  of  nest,  nests  with  the  probability 
of  paE  [0,  1],  In  this  way,  new  random  solutions  (new  nests)  can  replace  host  nests. 

Let's  discuss  the  rules  above,  (CS)  have  been  implemented  as  follows.  The  solution  is  the  whole  egg  in  the  nest. 
Thus,  each  cuckoo  can  place  only  one  egg  in  a  nest  in  its  original  form,  although  each  nest  can  contain  multiple  eggs 
representing  a  range  of  solutions,  in  general.  Therefore,  there  is  the  possibility  that  each  bird  and  hatches,  put  only  one  egg 
in  the  nest  in  the  same  original  form,  giving  us  that  each  nest  contains  eggs  of  multiple  qualities  representing  a  range  of 
solutions,  as  a  general  idea.  Mathematically,  these  types  of  trivial  problems  of  transformation  are  minimized  to  the 
maximum  of  problems  with  respect  to  the  micro-equation  (f  (x) )  =  max  (-f  (x)).  In  contrast  to  the  objective  function,  this 
function  is  now  referred  to  as  a  fitness  function  [8], 

Nature's  axioms,  animals  search  for  food  in  a  random  or  semi-random  fashion.  For  random  search,  the  animal 
search  path  is  randomized  without  any  discrimination  because  the  next  step  depends  on  (location  -  current  status  - 
probability  of  moving  to  your  next  location).  The  chosen  direction  implicitly  depends  on  the  probability  of  the 
mathematical  model  [30].  From  all  of  the  above,  the  purpose  of  the  CS  algorithm  can  be  written  clearly  is  to  look  for  x  that 
reduces  the  function  of  the  objective  function  f  (x).  Through  these  rules,  we  conclude  that  for  the  cuckoo  nest  i,  the  new 
generation  solution  can  be  defined  within  the  k-th  iteration  as  [30]. 

xf+1  =  xf  +  a  ®  Levy  (/?)  (1) 

To  illustrate  the  idea,  we  have  to  give  each  symbol  a  symbol: 

a  is  the  constant  step  of  the  method  of  Flying  Levy  possible  using  its  value  a=  1.5. 
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©  Product  code  which  means  double  entries. 

After  passing  through  the  research  on  the  behaviour  of  random  search  into  the  nature  of  animals  with  proof  of 
behaviour  characteristics  meet  Levy  flight.  We  can  equate  the  random  walk  and  the  amount  of  repetition  following  by  the 
Levy  flight  using  the  distribution: 

Levy  ~  u  —  k~ P  ,  0  <  (3  <  3  (2) 

The  variable  u  is  certainly  the  normal  distribution  and  the  number  of  repetitions  is  k.  How  many  times  do  we 
randomly  generate  these  numbers,  obey  Levy  flights,  and  calculating  the  amount  and  direction  of  the  step  in  a  random 
manner  using  the  Levy  distribution  given  in  the  two  research  papers  [6],  [18]. 

Which  has  (contrast  with  mean)  infinite.  Here,  successive  sequence  /  cuckoo  steps  are  essentially  the  random 
walking  process  that  corresponds  to  the  energy  distribution  of  the  step  length  law  with  a  heavy  tail.  We  can  say  that  if  the 
cuckoo's  egg  is  very  similar  to  host  eggs,  it  is  unlikely  to  detect  this  coco  egg,  and  therefore  fitness  must  be  linked  to 
different  solutions.  Therefore,  it  is  recommended  to  perform  a  random  path  in  a  biased  manner  with  some  random  step 
sizes  ([9],  [30]).  In  Figure  1,  in  the  Flowchart  mode  we  show  the  Cuckoo  search  path  to  access  the  function  formula: 


Figure  1:  Illustrated  the  Cuckoo  Search  Algorithm  [1] 


One  of  the  most  important  features  of  this  algorithm  is  its  simplicity.  If  we  compare  it  with  the  rest  of  the 
algorithms  that  are  categorized  by  metaheuristics  algorithms  based  on  population  or  proxy  (PSO)  and  (HS),  there  is 
essentially  only  one  parameter  pK  in  CS  (regardless  of  population  size  n).  Therefore,  very  easy  to  implement.  After 
speaking  in  general  about  the  cuckoo  algorithm  and  slow  performance  we  give  an  introduction  to  some  types  of  Barrier 
function  and  define  the  details  of  its  formulas  in  the  following  section: 
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2.2.  Barrier  Function 

The  use  of  the  SUMT  technique  has  recently  been  popular,  with  the  original  restricted  problem  becoming  a  series 
of  unrestricted  problems.  This  was  a  natural  approach  because  the  non-codified  reduction  techniques  were  well  developed. 
As  research  progressed  in  optimization  algorithms,  other  methods  were  shown  to  be  more  efficient  and  reliable  for 
"typical"  optimization  problems.  However,  the  larger  the  problem,  the  more  appropriate  methods  can  become  numerically 
ineffective.  Thus,  especially  for  structural  optimization  using  rounding  techniques,  the  new  appearance  of  SUMT  is 
appropriate  [28]. 

The  famous  method  in  community  optimization  metaheuristic  to  deal  with  constraints  is  to  use  the  Barrier 
function.  The  main  idea  of  this  method  is  to  convert  the  problem  of  restricted  improvement  into  an  unrestricted  problem  by 
adding  numerical  coefficients  for  both  the  objective  function  and  the  constraints  given  to  the  question  depending  on  the 
method  used  and  the  number  and  type  of  constraints  to  the  problem.  This  technique,  known  as  the  Barrier  method,  is  one  of 
the  famous  methods  of  dealing  with  evolutionary  algorithms.  The  current  work  has  been  using  the  same  method  [29]. 

2.3.  Modified  Frish's&  Carroll's  Barrier  Functions 

To  view  the  nonlinear  programming  problem  (NLP)  with  the  inequality  limitations,  we  give  the  following  sample: 
minimize  f0  (x) 

subject  to  q(x)  >  0,  i  =  1, ... ,  m  (3) 

Where  x  £  Rn,  C[(x):  Rn  — >  R,  are  real-valued  functions.  The  classical  (linear)  Lagrangian  of  (NLP)  is  defined  by: 

L(x, u)  =  f0(x)  -  Z^rViCj,  i=l, . ,m  (4) 

There  are  factors  that  are  important  components  in  explaining  the  optimum  conditions  for  NLP  and  the 
combination  of  algorithms  to  solve  NLP.  There  are  two  forms  of  nonlinear  functions  that  can  be  addressed  as  they  are: 

•  The  convex  programming,  through  the  writing  of  the  saddle-point  theorem  based  on  the  classical  Lagrange 
method  and  the  binary  methods  used,  the  solution  L(x,  vk)  of  the  mini-vk  can  be  updated. 

•  Nonconvex  programming,  L  (x,  vk)  is  often  the  most  important  form  of  a  non-convex  shape  for  vk,  which  is 
closer  to  v*  and  x  in  the  vicinity  of  x*,  where  x  (v*)  is  Kuhn-Tucker  (NLP),  which  creates  complexity  in 
numerical  issues. 

It  is  possible  to  say  that  those  who  solved  this  question,  the  researcher  Hestenes  through  his  research  [11]  and 
Powell  [23]  Lagrangian  extended  questions  of  addition  to  the  constraints  of  equality  and  Rockafellar  [24],  [25]  (Equality 
and  inequality),  but  the  convex  programming,  Polyak  and  Teboulle  [20]  discussed  a  different  class  of  Lagrange  based  on 
the  Log-Sigmoid  function.  Polyak  [22]  then  developed  the  NR  method  and  gave  the  approach  of  the  method  in  detail.  In 
[2]  they  presented  four  types  of  hybrid  structures  for  the  function  of  Polyak's  Log-Sigmoid  with  a  firefly  algorithm  at  a 
time  of  this  year.  As  for  nonconvex  programming,  a  class  of  nonlinear  langarians  is  usually  given  with  the  limitations  of 
inequality  to  the  question,  which  then  gives  us  the  unencumbered  saddle  points  by  the  Mangasarian  [15].  After  looking  at 
the  lower  pth  function  given  by  Charalambous  [5];  the  researcher  Bertsekas  [3]  modified  the  exponential  Lagrangian 
function  as  follows: 
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F(x,  u,  k)  =  f„(x)  +  [l  -  ekci«] 

Polyak  [21]  gave  two  modified  barrier  functions,  namely,  modified  Frish’s  function 


F(x,  u,k)  = 


f0(x)  -  ^2SiViln(l  +  roCi(x)),x  £  intSk 
+oo,  x  £  intSk 


(5) 


(6) 


(7) 


and  modified  Carroll’s  function: 

C(x  u  k)  =  [fo(x)_“^iVi(l-[l+(BCi(x)]-1),x£  intSk 
l  +oo,  x  £  int  Sk 

where  k  >  0  is  parameter  and  Sk  =  {x|  1  +  a>Ci(x)  >  0,  i  =  1, ... ,  m}. 

In  this  article,  we  tried  to  make  the  most  of  these  functions  in  (6)  and  (7).  We  proposed  a  fractional  function 
consisting  of  a  numerator  and  denominator  as  given  in  the  following  equation: 

u  Modified  Frish’s  function 


,  ,  =  New  Rational  Barrier  Function 

(l-u)Modified  Carrolrs  function 

Where  is  the  scalar  between  0<u<l. 

3.  MIXED  CS  AND  RATIONAL  BARRIER  FUNCTION  CRB 


(8) 


When  dealing  with  general  functions  (linear  -  nonlinear),  it  is  worth  noting  the  difficulties  in  the  question  when 
finding  the  best  solution  for  them  and  this  is  followed  by  formulas  and  research  in  finding  the  best  ways.  In  this  article,  we 
combined  two  important  types.  The  first  type  of  treatment  of  the  question  and  its  transformation  from  a  functionally 
restricted  constraint  and  nonlinear  constraints  by  the  function  of  the  barrier  with  the  egg  distribution  method  to  the  gums 
gave  a  good  and  appropriate  formula  for  solving  this  type  of  issue.  Take  a  look  at  Figure  2  we  notice  that  the  new 
algorithm  method  based  on  the  division  of  two  functions  developed  from  the  functions  of  the  barrier  (Frish’s  &  Carroll’s 
functions)  and  multiplied  by  a  coefficient  is  the  highest  level  at  1  and  the  least  possible  to  take  is  the  0  in  the  significance 
and  strength  of  the  fractional  formula.  Cuckoo  to  get  better  behaviour  and  approach  the  optimal  function  of  the  problem 
when  solving  non-linear  functions  multiple  constraints  and  complex  engineering  functions.: 


www.tjprc.ors 


editor@tjprc.  org 


18 


Rana  Z.  AI-Kawaz,  Hilda  I.  Ahmed  &  Abbas  Y.  Al-Bayati 


Figure  2:  The  New  Algorithm  CRB 


4.  STRUCTURAL  ENGINEERING  OPTIMIZATION 

In  this  section  we  will  show  the  efficiency  of  our  new  algorithm  when  compared  with  15  nonlinear  functions 
given  in  the  Appendix  within  the  article  and  a  number  of  complex  structural  optimization  problems,  and  sometimes  there 
are  no  optimal  solutions  for  attention.  In  order  to  find  out  how  the  CRB  algorithm  performs  against  the  CS  algorithm  about 
3  problems  solve  standard  structural  engineering  testing. 

Case  1:  Spring  Pressure 

The  geometric  questions  given  in  [4]  consist  of  the  following  form  and  data,  namely  weight  reduction  for  pressure 
pressure  (Figure  3)  controlled  by  minimum  deviation  limits,  shear  stress,  frequency  of  mutation,  outer  diameter  perimeter 
and  design  variables.  There  are  three  design  variables:  the  XI  wire  diameter,  the  average  diameter  of  the  X2  coil,  and  the 
number  of  active  X3  files: 
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Case  2:  Speed  Reducer 

The  speed  reduction  question  [4]  shown  in  the  Figure  4,  with  variables  defined  as  follows: 

After  the  face  X[,  the  tooth  unit  X2,  the  number  of  teeth  on  the  wing  X3,  the  length  of  the  column  1  of  the  X4 
axillary,  the  length  of  column  2  between  the  bearings  X5  and  diameter  (column  1)  X6  and  the  diameter  of  the  first  X7  (all 
variables  are  continuous  except  X3  which  is  integer).  The  weight  of  the  reducer  should  be  reduced  according  to  the  gears  of 
the  gear  teeth,  the  surface  stress,  the  occasional  deviations  of  the  columns  and  the  pressures  in  the  column. 


I 


r 


.V, 

T. 


Figure  4:  Speed  Reducer 

Case  3:  Welded  Beam 

The  last  engineering  questions  are  the  welded  package  at  the  lowest  cost  [4].  Figure  5  shows  the  beam  structure  of 
the  beam  A  and  the  weld  required  to  load  it  to  member  B.  The  question  is  to  find  the  minimum  manufacturing  cost  with  its 
four  components: 

Xi,  X2,  X3,  X4,  x,  bending  stress  in  the  beam  a,  buckling  load  on  the  bar  Pc,  and  end  deflection  on  the  beam  5. 


Figure  5:  Welded  Beam 

The  values  in  tables  starting  from  Table  1  to  Table  4  show  the  results  of  the  best  performance  of  the  new  algorithm  in 
terms  of  the  number  of  tools  described  in  the  tables  below: 
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Table  1:  Comparing  CS  Algorithm  against  CRB  in  2D  models  for  Number  of  nests=25 
(Total  Number  of  Iterations=10000,  Time=2000) 


Function 

CS 

CRB 

fmin 

Mean  Best  Nest 

fmin 

Mean  Best  Nest 

FI 

1 ,400000000000000e+0 1 0 

3.179203962030118 

1.380687544362068e+003 

1.792844905347926 

F2 

1 ,500000000000000e+0 1 0 

5.155094546754263 

1.000000000000000e+010 

3.643501189377146 

F3 

1 .000000000000000e+010 

2.104849965082111 

0.199946282778281 

1.599989645305540 

F5 

1 ,300000000000000e+0 1 0 

5.372292496848714 

1 .000000000000000e+0 1 0 

3.229702777984627 

F6 

6.78429146418581  le+002 

2.169880798217203 

1.31 6848978579592e+002 

1.016848978579592e+002 

F7 

1 .000000000000000e+010 

2.810772173507388 

0.372381895049205 

1.120287140654425 

F8 

1 .000000000000000e+010 

2.506881169340502 

0.876789754634672 

1.406630442441339 

F10 

1.209992957493290 

0.799984915177992 

0.209984423794242 

0.499989304529111 

FI  1 

1 ,300000000000000e+0 1 0 

1.095447189578642 

1.000000000000000e+010 

0.127130574560910 

F12 

1 .000000000000000e+010 

2. 475735361799080 

1 .000000000000000e+0 1 0 

2.070593645908076 

F13 

1.156261555600451 

1.937498420682251 

0.156256259390069 

0.937499998704024 

F14 

1 .000000000000000e+010 

0.671277061080557 

7. 121486997773 123e+002 

0.112499978892667 

F15 

3.229286417411150 

0.945858903815962 

0.262572203391223 

0.189812314433039 

Table  2:  Comparing  CS  Algorithm  against  CRB  in  3D  and  Engineering 
Models  for  Number  of  Nests=25 
(Total  Number  of  Iterations=10000,  Time=2000) 


Function 

CS 

CRB 

fmin 

Mean  Best  Nest 

fmin 

Mean  Best  Nest 

F4 

1 .200000000000000e+0 1 0 

2.957017834544342 

1 .000000000000000e+0 1 0 

0.985248894239120 

F9 

1 .900000000000000e+0 1 0 

3.992116826916045 

1.602811434499159 

3.766666666666667 

Spring 

1 .000000000000000e+010 

5.361519719525025 

0.801571652220139 

0.737997242107562 

Speed 

3.2000000000000e+010 

7.832525695537344 

2.3620388 17325554e+003 

6.185714285714285 

Beam 

1.724852313132572 

3.229642967361158s 

0.309403299807856 

0.100000000000000 

Table  3:  Comparing  CS  Algorithm  Against  CRB  in  2D  Models  for  Number  of  Nests=25 
(Total  Number  of  Iterations=5000,  Time=100) 


Function 

CS 

CRB 

fmin 

Mean  Best  Nest 

fmin 

MEAN  Best  Nest 

FI 

1 ,400000000000000e+0 1 0 

3.179203962030118 

1 ,232828996529320e+008 

1.222217898395526 

F2 

1 ,500000000000000e+0 1 0 

5.155094546754263 

1.141 124068742142e+008 

1.026826803301498 

F3 

1 . 000000000000000e+0 1 0 

2.104849965082111 

1 .000000000000000e+0 1 0 

0.906582056092352 

F5 

1 .300000000000000e+0 1 0 

5.372292496848714 

1 .000000000000000e+0 1 0 

4.577859246219106 

F6 

6.78429146418581 le+002 

2.169880798217203 

1 .000000000000000e+0 1 0 

1.898848112696566 

F7 

1 . 000000000000000e+0 1 0 

2.810772173507388 

0.83791 389 1472630e+009 

1.380030554991092 

F8 

1 . 000000000000000e+0 1 0 

2.506881169340502 

1 .000000000000000e+0 1 0 

2.322684997789185 

F10 

1.209992957493290 

0.799984915177992 

1.1 30633 10067287 le+009 

0.647046115409421 

FI  1 

1 ,300000000000000e+010 

1.095447189578642 

1 .000000000000000e+0 1 0 

1.167527944425209 

F12 

1 .000000000000000e+010 

2.  475735361799080 

1 .000000000000000e+0 1 0 

1.257638143994151 

F13 

1.156261555600451 

1.937498420682251 

1. 208696388 185930e+004 

0.910749881706651 

F14 

1 .000000000000000e+010 

0.671277061080557 

1 .000000000000000e+0 1 0 

0.431168318976144 

F15 

3.229286417411150 

0.945858903815962 

0.6571 18204978861e+009 

0.199051380299188 

Table  4:  Comparing  CS  Algorithm  Against  CRB  in  3D  and  Engineering 
Models  for  Number  of  Nests=25 
(Total  Number  of  Iterations=5000,  Time=100) 


Function 

CS 

CRB 

fmin 

Mean  Best  Nest 

fmin 

Mean  Best  Nest 

F4 

1 ,200000000000000e+0 1 0 

2.957017834544342 

1 .000000000000000e+0 1 0 

1.674673521443532 

F9 

1 . 900000000000000e+0 1 0 

3.992116826916045 

0.602811434499159 

3.766666666666667 

Spring 

1 .000000000000000e+010 

5.361519719525025 

0.393036466676572 

0.674020005115085 

Speed 

3 . 2000000000000e+0 1 0 

7.832525695537344 

2.3620388 17325554e+003 

6.185714285714285 
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Beam 

1.724852313132572 

3.229642967361158s 

0.309403299807856 

0.100000000000000 

5.  CONCLUSIONS 

The  results  of  nesting  and  minimum  function  values  in  the  Tables  1-4  showed  that  functions  (1-18)  have  a  more 
appropriate  reduction  than  the  use  of  the  Cuckoo  algorithm  alone  on  these  functions.  We  have  reviewed  through  the  tables 
that  when  reducing  frequency  algorithm  performance  is  better  than  when  increasing  the  number  of  repetitions  and  later  we 
proved  that  the  new  method  has  a  higher  capacity  than  the  algorithm  of  the  cuckoo  alone  and  its  strong  impact  on 
engineering  functions.  Calculations  were  done  using  the  Matlab  2011  program  and  with  a  calculator  Core  i5. 

6.  APPENDIX 

1- rnin  f(x)  =  (xx  —  3)2  +  (x2  —  2)2 

s.t  xf  —  x|  +  5  =  0 
x1  +  2x2  —  4  <  0 

2- minf(x)  =  x2  —  XjX2  +  x| 

s.t  x2  +  x|  -  4  =  0 
2x1  —  x2  +  2  <  0 

3- minf(x)  =  (xx  -  2)2  +  (x2  -  l)2 

s.t  x-l  —  2x2  =  —1 
—  —  +  x|  +  1  >0 

4- minf(x)  =  x2  +  2x|  x3  +  2x3 

s.t  xf  +  x2  +  x|  =  4 
xi  —  x2  +  2x3  <  2 

5- minf(x)  =  e1  —  xxx2  +  x| 

s.t  x\  +  x|  =  4 
2x1  +  x2  <  2 

6- minf(x)  =  x2  —  3x1x2  +  4 

s.t  —  2x1  +  x|  =  5 
5xx  +  2x2  >  18 

7- minf(x)  =  —  e~Xl~X2 

s.t  x2  +  x2  —  4 
xt  —  1  >  0 

8- minf(x)  =  -x2  +  2xxx2  +  x2  -  e-Xl_X2 

s.t  x2  +  x2  —  4  —  0 
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Xt  +  X2  <  1 

9- minf(x)  =  — x1x2x3 

s.t  20  —  x1  >0 
11  —  x2  >0 
42  —  x3  >  0 

72  —  x1  —  2x2  —  2x3  >  0 

10- minf(x)  =  (xx  —  l)2  +  x2  —  2 

s.t  x2  —  xx  =  1 
x1+x2  >  2 

11- minf(x)  =  xx2  +  x22 

s.t 

.  x1  —  3  =  0 
x2  —  2  <  0 

12- minf(x)  =  ^(x^  +  x22)  -  cosCxJ  cos  (^|)  +  1 

s.t  xt  —  3  =  0 
x2  —  2  <  0 

13- minf(x)  =  (xx  -  2)2  +  ^x22 

s.t  2xx  +  3x2  =  4 
Xi  -  ^x2  <  1 

14- minf(x)  =  — xxx2 

s.t  20xx  +  15x2  —  30  =  0 
— +  x22  -  1  <  0 

4  Z 

15- minf(x)  =  xx4  —  2xx2x2  +  xx2  +  xxx22  —  2x1  +  4 

s.tXj2  +  x22  -  2  =  0 
0.25  xx2  +  0.75x22  -  1  <  0 
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